What can we learn about neutron stars from gravity-wave observations?* 
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In the next few years, the first detections of gravity-wave signals using Earth-based interferomet- 
ric detectors will begin to provide precious new information about the structure and dynamics of 
compact bodies such as neutron stars. The intrinsic weakness of gravity-wave signals requires a 
proactive approach to modeling the prospective sources and anticipating the shape of the signals 
that we seek to detect. Full-blown 3-D numerical simulations of the sources are playing and will play 
an important role in planning the gravity-wave data-analysis effort. I review some recent analytical 
and numerical work on neutron stars as sources of gravity waves. 
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I. INTRODUCTION 

In the course of the next decade, the inception of 
gravity-wave (GW) astronomy will open an exciting 
new window on the physics of compact, strongly grav- 
itating objects such as neutron stars (NSs) and black 
holes (BHs), providing information complementary to 
that available from electromagnetic and neutrino obser- 
vations, and plausibly producing important insights into 
unsolved questions such as the equation of state (EOS) of 
matter at nuclear densities, and the mechanisms behind 
gamma-ray bursts. Detailed reviews of GW sources, of 
the expected event rates, and of the physics thai, these 
sources could teach us are available elsewhere, aaB but I 
shall list briefly the most promising astrophysical systems 
from which we could learn about NSs using Earth-based 
GW interferometers such as LIGO and VIRGO. 

1. NS-NS and NS-BH binaries in the last few min- 
utes of their inspirais. For a long time, these in- 
spiraling systems have been the prototype for the 
category of the short-lived chirp signals detectable 
using Earth-based interferometers. The reason, of 
course, is that NS-NS binaries have actually been 
observed in our galaxyjj but also that the part of 
the inspiral accessible to the interferometers (with 
GW frequencies between 40 and 1000 Hz) sits well 
before the final merger of the binary, so it is de- 
scribed very accurately by the well-developed post- 
Newtonian equations for point masses. a The suc- 
cessful observation of GWs from these events will 
teach us about the masses, spins and locations of 
NSs, but not about their internal structure. 

On the contrary, the detection of GW from the end- 
point of NS-BH inspirais should produce detailed 
information about NS structure and EOS. For a 
wide range of binary parameters, the NS will be 
torn apart by the tidal field of the BH well be- 
fore the final plunge into the hole, and the tidal- 
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disruption waves will be well inside the frequency 
range of good interferometer sensitivity.El NS-BH 
binaries have also-, been proposed as engines for 
gamma-ray burstso and as suitable environments 
for the production of heavy nuclei in r-processes.u 
I will discuss these systems more extensively in 
Sec. |n] of this paper. 

2. Rapidly spinning, deformed NSs. This class in- 
cludes the known and unknown radio pulsars (when 
their gravitational ellipticity is high enough to pro- 
vide strong GWs), and the systems known as low- 
mass X-ray binaries (LMXB), where the NS is ac- 
creting matter and angular momentum from a com- 
panion, but instead of increasing its rotation, it is 
locked into spin periods of about 3 ms; it is conjec- 
tured that the angular mome,atuni being accreted 
is lost to the emission of GW.Ia 

To detect rapidly spinning NSs, it will be neces- 
sary to integrate the GW signal for times up to 
several months, so the Doppler frequency modula- 
tion caused by the movement of the Earth around 
the sun will make it mueh harder to detect pre- 
viously unknown sources.li2l At the same time, the 
shape of this modulation will make it possible to 
obtain the position of the source in the sky,E2l and 
to match the GW source with one of the objects 
known from electromagnetic observations. 

If any GWs are detected, their features would be 
very informative, in particular when examined in 
correlation with electromagnetic signals from the 
same source. For instance, the ratio of the GW 
frequency to the NS angular frequency could iden- 
tify the nature of the inhomogeneities that give rise 
to the GW emission, and the evolution of the GW 
amplitude and frequency could provide interesting 
data about NS physics such as crust structure and 
dynamics, crust-core interactions™ magnetic fields, 
viscosity, superfluidity, and more.cl 

3. Proto-neutron stars. Finally, NSs could be ob- 
served as the rapidly spinning, strongly asymmetric 
remnants of stellar-core collapse, or as the proto- 
NSs produced by the accretion-induced collapse of 
white dwarves. Proto-NSs that spin very fast can 
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hang up centrifugally at a stage where their ra- 
dius is still large compared to that of the final NS. 
Such a configuration might be unstable to a bar 
mode, giving rise to an elongated object that would 
emit very strong GWs.ta The newborn NSs might 
also develop a GW-induced instability in their r- 
modes.H3ll3 I wi ll d iscuss this possibility more ex- 
tensively in Sec. 
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For all these systems, GWs would provide informa- 
tion complementary to that made available by neu- 
trino observations, focusing on the density struc- 
ture and asymmetry of the collapsing core rather 
than on its thermal structure. 

While we wait eagerly for the first detections, numeri- 
cal simulations of GW sources are beginning to provide 
precious insight into the effect of NS structure and dy- 
namics on the GW signals. The purpose of this paper 
is to explore how simulations are helping and directing 
the theory and practice of GW data analysis. To do so, 
I will single out two promising GW sources (NS tidal 
disruption in NS-BH binaries, in Sec. H NS r-modes, m 
Sec. [II) that are currently being attacked with numerical 
simulations, and on which I have first-hand experience. 



II. NEUTRON-STAR TIDAL DISRUPTION AS 
A PROBE INTO THE EQUATION OF STATE OF 
DENSE NUCLEAR MATTER 

With event rates between 10 -7 and 10 -4 per year in 
our galaxy, NS-BH mergers are one of the standard GW 
sources for second-generation interferometers (LIGO-II 
should be able to detect these systems out to 650 Mpc, 
yielding 1-1500 observations per yearn). The waveforms 
generated by these events will contain two kinds of in- 
formation. The early part of the inspiral (during which 
the NS and BH are still relatively distant, and the dy- 
namics can be described accurately by post-Newtonian 
equations of motion in the point-mass approximation) 
will tell us about the masses and the spins of the NS and 
BH. The late part of the inspiral, depending on the bi- 
nary parameters, can see the BH tidal field become so 
strong that it disrupts the NS on a dynamical timescale; 
physical intuition then suggests that the details of the dis- 
ruption process, as encoded in GWs, should carry useful 
information about the interna] structure of the NS, and 
in particular about its EOS.E3 



A. A simple analytical model for NS tidal 
disruption 

The prospects for extracting this information from the 
GWs that could be measured from a realistic event have 
been evaluated by the present author|j using a crude 
model that however accounts for most of the relevant 
physics. The simplest possible representation of a NS 
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FIG. 1: GW frequency at tidal disruption, /td(m, 7?, M, a), as 
a function of NS radius R. Here we have set the NS mass to 
1.4 Mq , we have indicated the BH mass on the curves, and we 
have omitted the dependence on the BH spin a, which gives 
negligible corrections. The BH spin parameter a is neverthe- 
less important to decide whether the NS will actually disrupt 
before it begins its rapid plunge into the BH. Estimating the 
onset of the plunge as the naive 6M ISCO, we find that to 
have disruption we need a/M to assume at least the value in- 
dicated by the big dots in this figure. Negative a/M denotes 
counterrotating orbits. 



inspiraling into a BH is a quasi-equilibrium sequence of 
relativistic Roche-Riemann ellipsoids. These ellipsoids 
are equilibrium configurations of a self-gravitating, poly- 
tropic, Newtonian fluid, moving on circular, equatorial 
geodesies in the Kert-spacetime, and subject to the BH 
relativistic tidal field.E3 For these configurations, once the 
binary parameters m, M , a and r (respectively, the NS 
and BH masses, the BH spin, and the binary separation) 
are fixed, there is only one free parameter, corresponding 
to the NS radius R. So we will take R as a representative 
of the uncertainty in the NS EOS. 

As it inspirals toward the BH, a NS with parameters m 
and R would be represented by the appropriate Roche- 
Riemann ellipsoid at each separation r, until we reach 
a critical r cr , beyond which no more equilibrium config- 
urations exist. We identify this end of the equilibrium 
sequence with the onset of dynamical tidal disruption, 
and from r cr we obtain the GW frequency at tidal dis- 
ruption, /td = /td(w, i?, M, a). This function is shown in 
Fig. for the standard NS mass m = 1.4 Mq, and for a 
variety of BH masses. Figure |l| should be read as follows: 
given a NS-BH signal, choose the curve corresponding to 
the BH mass (as estimated from the frequency evolution 
of the early inspiral signal); on the horizontal axis, locate 
the frequency at which tidal disruption begins (as esti- 
mated from the late part of the GW signal); then read 
off the estimated NS radius on the vertical axis. 

Figure [l] suggests that the GW frequency at tidal dis- 
ruption depends strongly on the NS radius, and that the 
disruption waveforms lie in the band of good interferon!- 
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eter sensitivity for the advanced interferometers such as 
LIGO-II. It follows that, in principle, we could use the 
waveforms from a NS tidal-disruption event to measure 
both the NS mass and the NS radius. However, we still 
need to know just how well we could measure them. 



B. The matched-filtering paradigm for 
gravity-wave data analysis 

To answer this question, we need to invoke the theory 
of matched-filtering parameter estimationSJl The general 
idea is that GWs will be detected by correlating the mea- 
sured signal |s) to a bank of theoretical templates {\hi)}, 
which represent our best approximation of the realistic 
GW signal, as it could happen for a variety of binary 
parameters. If the match (s\hi) (the correlation between 
\hi) and \s)) is much higher than the match (hi\n) that 
the template would give, on the average, with noise alone, 
then we claim that we have a detection. To know how 
well we can estimate R, we ask how probable it is that 
a particular realization of detector noise would lead us 
to mistake the template \hu) with the nearby template 
the answer is given in terms of the match 
(hn\h R+ AR)- 

To compute this match-we construct a bank of signal 
templates that differ onlyc3 in the GW frequency at the 
onset of tidal disruption. Because we are not interested 
in modeling accurately the relativistic dynamics of the 
binary, but only the effects of tidal disruption, we gen- 
erate our waveforms-Jjrom simple quadrupole-governed 
Newtonian inscirals,E3 cutting off the signal, more or 
less abruptly,c3 when the instantaneous GW frequency 
reaches ftd- Computing the match between nearby tem- 
plates, we estimate the granularity to which / t( j can be 
measured for a given signal strength (inversely propor- 
tional to distance); we can then use our equation for f(R) 
to propagate the errors to the NS radius. The final result 
of this exercise is that, employing, advanced Earth-based 
interferometers such as LIGO-II0, we should be able to 
measure the NS radius to 15%, for tidal-disruption.events 
at distances that yield about one event per yearfl This 
estimated error should be compared to the error (about 
a factor of two) in the measurements of the NS thermal 
radius. 



C. Inverting the mass— radius relation 

What can we do once we get a few points on the m(R) 
curve? We can try to solve for the EOS of dense nu- 
clear matter in NSs. The relativistic model of nonro- 
tating NSs may be considered as a mapping from the NS 
EOS, through the Oppenheimer-Volkoff (OV) equations, 

— = Anr 2 p d(> = (pi p) m + 4?rr3p (1) 
dr dr r(r — 2m) 



to macroscopic NS quantities such as mass and radius. 
Let us work through this mapping in the case of a poly- 
tropic EOS, p — p(p). First, we set the central density 
p c : then, we solve the OV equations and compute the 
NS radius R = R{p c ) and mass m = m(p c ); finally, we 
eliminate p c from these two equations, completing our 
mapping of the EOS p{p) into the mass-radius relation 
m(R). 

LindblomEHl has shown how to invert the OV mapping 
using even a few m(R) data points. Suppose that p(p) 
is known up to a certain density p max from other obser- 
vations and experiments. Now start with the least dense 
NS, and integrate the OV equations backward, from the 
surface of the star [where we know R and m(R)] down 
to the radius where p — p ma x- We are then left with a 
stellar core of known mass and radius, and we can use an 
analytic approximation for the solution of the OV equa- 
tions to get p c and p c ; we add the point p c (Pc) to the 
EOS, and repeat for the next NS. The result is p(p) as 
a piecewise linear law. This analysis can be generalized 
to include rotationally deformed models, and to account 
for the statistical uncertainty in the m(R) data. Re- 
cently, HaradaE!l has shown how other macroscopic pa- 
rameters of NSs (such as moments of inertia, baryonic 
masses, binding energies, gravitational redshifts) can be 
used in this same framework to recover information about 
the EOS. 

Finally, Saijo and Nakamurac3 have suggested that 
it might be possible to measure the NS radius directly 
from the spectrum of the GWs emitted in NS-BH coa- 
lescences. These authors have used BH perturbation the- 
ory to compute the spectrum of GW emitted by a disk 
of dust inspiraling into a rotating BH. When the radius 
i?disk is larger than the wavelength of the quasi-normal 
modes of the BH, the spectrum acquires several peaks 
with separation oc -R^isk> irrespective of M and a. Saijo 
and Nakamura conjecture that the same structure would 
be visible in the spectrum of GW signals from NS-BH 
binaries, providing direct information about the radius. 
However, several issues are left unaddressed: in particu- 
lar, the particles of the disk move along geodesies, while 
the fluid of a NS would be strongly constrained by the 
gravitation and pressure of the star (except, perhaps, in 
the regime of severe tidal disruption); furthermore, for 
coalescence events that happen at realistic distances the 
signal strength might be too low to let us resolve the 
form-factor structure in the spectrum. 



D. Numerical simulations of neutron-star tidal 
disruption 

Summarizing, the preliminary analysis carried out by 
the present author suggests that NS tidal-disruption 
events have much to teach us about the NS EOS. Jjhis 
prediction is confirmed by Newtonian simulationsE3 of 
NS-BH systems, carried out using both smooth parti- 
cle hydrodynamicsocj and Eulerian techniquesQCJ, which 
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show that the ultimate fate of the system (complete or in- 
complete disruption; presence of accretion rings or tidal 
tails; features of the GW shutoff) depends strongly on 
the stiffness of the EOS. However, detailed relativistic 
numerical simulations are still needed to confirm these 
prospects, and will be essential as a foundation to inter- 
pret any tidal-disruption waveforms that might be mea- 
sured in reality. 



III. GRAVITY WAVES FROM THE r-MODES 
OF YOUNG NEUTRON STARS 

All rotating stars possess a class of circulation modes 
(r-modes) that are driven toward instability by gravita- 
tional radiation reaction; in hot, rapidly rotating young 
NSs, this destabilizing effect might be so strong that it 
dominates viscous dissipation. Once an r-mode achieves 
sufficient amplitude, the star is quickly spun down as-ap. 
gular momentum is lost to gravitational radiation. EjEia 
For stars with initial angular velocity £1 ~ 1000 Hz, the 
timescale for r-mode growth is ~ 40 s. In recent years, r- 
modes have attracted considerable interest as a possible 
explanation for the failure to observe very rapidly spin- 
ning pulsars, and as a promising GW source for Earth- 
based detectors. 

However, the astrophysical relevance of r-modes is still 
in doubt, pending judgment on two separate issues. First, 
the instability of r-modes (discovered by analyzing the 
linearized Euler equations for perfect fluids) might not be 
confirmed after all the complicated physics that occurs 
in NSs is taken into account, including relativistic effects, 
physical-EOS effects, solid crust effects, magnetic fields, 
rotation laws, and so on. Second, the amount of angular 
momentum removed from the star and the strength of the 
GW radiation emitted depend critically on the maximum 
amplitude that can be reached by the r-mode; but the 
growth of the r-mode might be limited to a very small 
saturation amplitude by the effects of magnetic fieldscJ, 
or by the leakage of energy to other (damped) modes by 
way of nonlinear hydrodynamical couplings. 

There have been several attempts to investigate 
the nonlinear dynamics of the r-modes, by -means 
of second-order Lagrangian perturbation theory,E3 and 
of relativistic pUumerical simulations in the Cowling 
approximation.^ In addition, Lindblom, Tohlinc and the 
present author have carried out Newtonian numerical 
simulations that include radiation reaction as an effec- 
tive force. We now briefly review this wor|k-,rbut we invite 
the reader to refer to the original articlesE3c3 for further 
details. 



A. Numerical evolutions of nonlinear r-modes in 
neutron stars 

We solved the Newtonian Euler and Poisson equa- 
tions on a 128 x 64 x 128 cylindrical grid, using a 2nd- 
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FIG. 2: Evolution of the dimensionless r-mode amplitude a. 
Amplitudes of order one imply mode velocities comparable to 
the local velocity of rigid rotation. Time is measured in units 
of the initial rotation period Po of the NS. 



order-accurate finite-difference code developed at LS 




to tackle a variety of astrophysical problems. The LSU 
code was parallelizedtEl in the domain-decomposition 
paradigm using the well known MPItS library. We ran 
it on 16 nodes of the HP V-2500 supercomputers sX Cal- 
tech's Center for Advanced Computing Research.^ 

Our initial equilibrium configuration was a simple 
polytrope (with n = 1, M = 1.4 Mq, R = 12.5 km) ob- 
tained by solvuag selfconsistently the Bernoulli and Pois- 
son equations.^ We then augmented the rigid-rotation 
velocity field with a small-amplitude^-slow-rotation ap- 
proximation to an I = m = 2 r-mode,E3 

5v = a m (^ 2 Y£e^\ (2) 

where R and f2o are the radius and angular velocity of the 
unperturbed star, ao is the dimensionless r-mode ampli- 
tude, and is a vector spherical harmonic of the mag- 
netic type. Gravitational radiation reaction was added as 
an effective Newtonian forceLlI proportional to the deriva- 
tives of the current multipole moments of the star. 

Our computing budget set tight limits on the length 
of the simulation: one full rotation period of a rapidly 
rotating NS took up to ten thousand CPU hours. Conse- 
quently, we could not afford to follow the growth of the 
r-mode through one or more radiation timescales. In- 
stead, we increased the strength of the effective radiation- 
reaction force by a factor of about 4500, to bring down 
the r-mode growth timescale to values comparable to 
the rotation period. Even with this trick, we expect 
that the physical behavior observed in the simulations 
should still be realistic, because the nonlinear dynamics 
of the r-mode (including its couplings to other modes) 
should happen at the hydrodynamical timescale, which 
in our simulations is still much shorter than the radiation- 
reaction timescale. 

While computing the evolution of the star, we moni- 
tored the r-mode amplitude a by projecting the velocity 
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FIG. 3: Cresting waves at the surface of the NS, after the r- 
mode amplitude has reached its maximum. The waves create 
strong shocks that dump mode energy into thermal energy. 



of the fluid on Yfi , 



2|J 2 



(Q/R) J pr 6 dr' 



(3) 



where J 21 is the I = m = 2 current multipolc moment; 
similarly, we read off the mode frequency to as 
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dt 



(4) 



At the beginning of our evolution, these diagnostics re- 
covered the theoretical values ao and ujq within the ex- 
pected error. 

Figure || shows the evolution of the r-mode amplitude. 
For this simulation, the theoretical predictiont3 for the 
r-modc growth time was 10 initial rotation periods. The 
mode grew exponentially (in good accord with theory) 
until a ~ 2; then the mode started to be limited by some 
nonlinear process, the amplitude peaked at a = 3.35, and 
finally fell down very rapidly. A_movie of the crashing 
r-mode can be found on the webEB 

What nonlinear process was responsible for limiting 
the growth of the r-mode, and for causing its rapid 
demise? We discovered a clue when we examined the 
evolution of the energy and angular momentum in the 
course of the simulation. Even after the emission of an- 
gular momentum into GWs fell to zero, the star contin- 
ued to lose energy; something other than GWs must be 
responsible for this loss. 

We believe that we have found the culprit. To first or- 
der in the amplitude, the r-mode is only a velocity mode; 
to second order, however, there is an associated density 
perturbation, proportional to Y32, which appears as a 
wave with four crests (two in each hemisphere) on the 
surface of the star. As the amplitude reaches its maxi- 
mum, these propagating crests turn into large, breaking 
waves (see Fig. ||), and the edges of the waves develop 
strong shocks that dump kinetic energy into thermal en- 
ergy, killing the r-mode. 



B. Numerical simulations as a laboratory of stellar 
physics 

Our code provided a nice laboratory to perform several 
more evolutions and tests. 



1. We performed basic tests on the robustness of the 
code and of our diagnostics, evolving stars with dif- 
ferent angular velocities, with or without r-mode 
perturbations, and at different grid resolutions. 



2. We investigated the dependence of the satura- 
tion amplitude on the artificial amplification of 
radiation reaction. Our computing budget made 
it impossible to increase the radiation-reaction 
timescale; instead, we reduced it even more, finding 
that the r-mode saturates faster, but at essentially 
the same amplitude. 



3. We studied the unforced evolution of unit- 
amplitude r-modes, finding that they are essen- 
tially stable for as long as we could evolve them. 
These results are compatible with the relativistic 
evolution^Xalso unforced) performed by Stergioulas 
and Font £3 



A new picture for gravity-wave signals from 
r-modes 



In the. traditional scenario for GW emission from r- 
modesj^j the unstable r-mode would grow until it reached 
a dimensionless amplitude of about one, and then it 
would saturate and persist at that amplitude (for several 
months) until it would have lost most of its angular mo- 
mentum; during that time, the frequency of the r-mode 
would decrease in proportion with the NS spin. The 
prospective picture of r-mode GW signals that emerges 
from our evolutions is quite different. Most interesting, 
the r-mode spindown episodes are faster (only a few min- 
utes), and the GW frequency remains remarkably con- 
stant as the angular velocity of the star decreases. As 
a result, the search for r-mode signals changes from a 
pulsar-like search (which must account for the Dopplcr 
shifts generated by the movement of the Earth around 
the sun) to an easier chiip-like search, with encouraging 
prospects for detection.Cj 

Although recent studies of the effect ofj-opagnetic 
fields ,E3 and of exotic forms of bulk viscositycil suggest 
that the astrophysical significance of r-modes might be 
limited, the considerable uncertainty about the macro- 
scopic and microscopic structure and dynamics of NSs 
makes it reasonable to devote some resources to GW 
searches for r-mode signals similar to those predicted by 
the simulations described above. 
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